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High resolution eigenvalue spectra of several two- and three-dimensional superconducting mi- 
crowave cavities have been measured in the frequency range below 20 GHz and analyzed using a 
statistical measure which is given by the distribution of the normalized mode fluctuations. For 
chaotic systems the limit distribution is conjectured to show a universal Gaussian, whereas inte- 
grable systems should exhibit a non-Gaussian limit distribution. For the investigated Bunimovich 
stadium and the 3D-Sinai billiard we find that the distribution is in good agreement with this pre- 
diction. We study members of the family of limagon billiards, having mixed dynamics. It turns 
out that in this case the number of approximately 1000 eigenvalues for each billiard does not allow 
to observe significant deviations from a Gaussian, whereas an also measured circular billiard with 
regular dynamics shows the expected difference from a Gaussian. 
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I. INTRODUCTION 

One of the main research lines in quantum chaos is 
to investigate the statistics of energy levels of quantum 
systems whose classical counterpart is chaotic. A very 
popular class of systems are Euclidean billiards which 
are classically given by the free motion inside a domain 
Jl C R 2 with elastic reflections at the boundary dfl. The 
corresponding quantum billiard is described by the sta- 
tionary Schrodinger equation {% = 2m = 1) 



A^„(q) +k 2 n V>„(q) 







for q E SI 



(1) 



with Dirichlet boundary conditions VVi(q) = for 
qedQ. 

It has been conjectured that the energy level statistics 
of integrable systems can be described by a Poissonian 
random process [IJ], whereas classically strongly chaotic 
systems should obey the statistics of random matrix en- 
sembles like the GOE or GUE |^-^] . This implies for ex- 
ample that the nearest neighbor level spacing distribution 
is expected to show level repulsion for chaotic systems, 
in contrast to integrable systems which are expected to 
show level attraction. Surprisingly, this means that the 
statistics of classically chaotic systems are much more 
rigid than those of integrable systems. 

These conjectures have been tested successfully for sev- 
eral systems. However, there are exceptions for both 
integrable and chaotic systems. For example the so- 
called arithmetic systems, which are strongly chaotic, are 
found to have a level spacing distribution showing level- 
attraction similar to the Poissonian distribution [p|-^0[. 
Therefore a different statistics, the distribution of the 
normalized mode fluctuations has been proposed [pT|Jl2[ 



as a possible signature of quantum chaos. This statis- 
tics was first investigated in the integrable case for the 
eigenvalues of the Laplacian on a torus Jl3| and later in 
|14| the unnormalized fluctuations, possessing no limit 
distribution, have been studied for regular and chaotic 
billiards. For chaotic systems the limit distribution of 
the normalized mode fluctuations is conjectured |Tl 12 
to show a universal Gaussian, whereas integrable systems 
should exhibit a non-Gaussian limit distribution. This 
conjecture was tested successfully for several regular and 
chaotic billiard systems in [|l2|,| 15 1(| . 

By using two-dimensional microwave cavities quantum 
billiards can be simulated experimentally This 
is possible because of the equivalence of the stationary 
Schrodinger equation for quantum billiards and the cor- 
responding Hclmholtz equation for electromagnetic res- 
onators in two dimensions. In three dimensions the elec- 
tromagnetic Helmholtz equation is vectorial and cannot 
be reduced to an effective scalar form. Thus, it is struc- 
turally different from the scalar Schrodinger equation. 
Nevertheless the applicability of the statistical concepts 
developed in the theory of quantum chaos and random 
matrix theory is also given for such three-dimensional 
20|-22]. Therefore experiments with supercon- 



systems 



ducting microwave resonators provide in general eigen- 
value spectra of very high resolution for which an analysis 
of the distribution of the normalized mode fluctuations 
is interesting. 

The paper is organized as follows. In Sec. || the mode 
fluctuation distribution is introduced. The experimental 
set-up and the measurement of the eigenfrequencies us- 
ing superconducting microwave resonators are described 
in Sec. [II. In Sec. IV the analysis of the mode fluctu- 



ation distribution using the experimental data is carried 
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out. 



II. MODE FLUCTUATION DISTRIBUTION 

The analysis of the eigenvalue spectrum starts with the 
spectral staircase function 



N(k) = #{n | k n < k}, 



(2) 



which counts the number of energy levels below a given 
energy k. The mean behavior of N(k) is given by the gen- 
eralized Weyl's law [^3| , which reads for two-dimensional 
billiards with Dirichlet boundary conditions 



N(k) = ^k 2 -^k+C , 

4"7T 47T 



(3) 



where A is the area of the billiard, C the length of the 
boundary, and C takes curvature and corner contribu- 
tions into account. For three-dimensional electrodynam- 
ical billiards we have 



(4) 
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const 
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where V is the volume of the billiard, TZ the mean radius 
of the curvature over the surface a and 10 is the dihedral 
angle along the edges a [§4],|5| . 

In order to obtain a spectrum which is independent 
of the system specific constants, one considers the un- 
folded spectrum {x n := N(k n )} pq , p7| . Consequently 
the unfolded energy spectrum has a mean level spacing 
of unity. The counting function for the unfolded spec- 
trum will be denoted for simplicity again by N(x). Thus 
the fluctuating part of the spectral staircase function is 
given by 



AW(x) := N(x) - N(x) = N(x) - x 



(5) 



In the following we will assume that all spectra have been 
unfolded and that the systems are completely desym- 
metrized. 

The normalized mode fluctuations are given by 



W(x) 



N(i uc {x) 



(6) 



where D(x) is the variance 
S(c) 



D(x) := 



(c-l)x 



[iVfluc(y)] 2 dy , 



(7) 



with c > 1, and 3(c) is a correction necessary for inte- 
grable systems to obtain for W(x) a variance of one, see 
pa for details. The conjecture put forward in [pl]Jl2| 



can be formulated as follows (see [T^]): For bound con- 
servative and scaling systems the quantity W(x), Eq. 
possesses a limit distribution for x —* 00. This distribu- 
tion is absolutely continuous with respect to the Lebesgue 
measure on the real line, with a density P(W) defined by 



cT 



lim -; — - 

T^oo ( C — 1) T 



g{W{x))p{x/T) dx 



g{W)P{W) dW 



(8) 



where g(x) is a bounded continuous function, and p(t) > 
is a continuous density on [1, c] with — f, p(t) dt = 1. 

Furthermore, the limit distribution has zero mean and 
unit variance, 



00 00 
J WP(W)dW = , J W 2 P(W) dW = 1 



(9) 



Now, the basic conjecture of [|ll|jl2| reads as follows: 
If the corresponding classical system is strongly chaotic, 
having only isolated and unstable periodic orbits, then 
P(W) is universally a Gaussian, 



P(W) 



exp 



—W' 
2 



(10) 



In contrast, a classically integrable system leads to a non- 
Gaussian density P(W). 

For large classes of integrable systems it has been 
proven that the limit distribution is not Gaussian, see 
p6| , p8| -]3C|] and references therein. 

For chaotic systems this conjecture has been tested 
numerically in ]l^,^,^5), and experimentally in fl3l[ . A 
review and a detailed comparison between chaotic and 
non-chaotic systems is given in p6[ ]. A Gaussian dis- 
tribution for chaotic systems corresponds to maximum 
randomness for the fluctuations, whereas the mode fluc- 
tuation distribution for integrable systems is less random 

& ™ 

Berry's semiclassical analysis p2J gives results for the 

asymptotic behavior of the saturation behavior of the 
spectral rigidity. In a similar way one can determine the 
asymptotic behavior of D{x). For generic integrable bil- 
liards one has 



D(x) ~ const.\fx 



(11) 



For generic classically chaotic systems with anti-unitary 
symmetry (e.g. time-reversal symmetry) one obtains 



D(x) 



1 

2^ 



In a; + const. 



(12) 



As discussed in [[16[ one can extend the conjecture to 
such chaotic systems, for which iVfl uc (fc) is modulated by 
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a long range oscillation N\ ong (k). In this case one has 
to include the additional term in the unfolding process, 
{x n := N(k n ) + N\ ong (k n )}. This procedure has been 
used for instance in the case of the truncated hyperbola 
billiard, where a prominent contribution to N(k) is given 
by families of closed non-periodic orbits running into a 
boundary point where the curvature is discontinuous fl33| ; 
see |l2| , ^6| for the result of P(W) for this system. 

The same is also necessary for the stadium billiard, 
where a family of marginally stable orbits, the bouncing 
ball orbits (bbo) gives rise to a strong modulation Jl9],[34| . 
Taking the contribution of the bbo into account one ob- 
serves excellent agreement of P(W) with the Gaussian 
normal distribution l|16| . 



III. EXPERIMENT 

For a precise test of the distribution of the normalized 
mode fluctuations an accurate measurement of the res- 
onances of all investigated microwave billiards is neces- 
sary. Since 1991, we have experimentally studied several 
two- and three-dimensional systems using superconduct- 
ing microwave resonators of niobium. In Fig. 1 the shapes 
of some measured billiards with their dimensions are 
shown. Altogether five desymmetrized two-dimensional 
systems are investigated: a 7 = 1.8 Bunimovich stadium 
billiard ||^,[l9|,|36| , a circular billiard (not desymmetrized) 
and three members of the family of limagon billiards, 
which have first been theoretically studied as billiards in 
p7f . Their boundary is defined as the quadratic con- 
formal map of the unit disc onto the complex w-plane: 
w = z + Xz 2 , where A S [0,1/2] controls the chaotic- 
ity of the system. The billiards of the limagon family 
are also called Pascalian snails, their shape has already 
been mentioned by A. Durer in 1525 |38|]. We also an- 
alyzed a desymmetrized three-dimensional Sinai billiard 



The measurements were carried out in a LHe-bath 
cryostat. The billiards were excited up to a frequency 
of 20 GHz, the upper limit of the used HP8510B vector 
network analyzer, using four capacitively coupling dipole 
antennas sitting in small holes on the niobium surface 
and penetrating up to a maximum of 1 mm into the 
cavity to avoid disturbances of the electromagnetic field 
inside the resonators. Using one antenna for the excita- 
tion and either another or the same one for the detec- 
tion of the microwave signal, we were able to measure 
the transmission or the reflection spectrum of the res- 
onators. The spectra were taken in 10 kHz-steps and 
the measured resonances have quality factors of up to 
Q Ri 10 7 and signal-to-noise ratios of up to S/N ~ 70 dB 
which made it easy to separate the resonances from each 
other and from the background even in the higher fre- 
quency range, where the level density strongly increases. 
Especially in the case of the circular billiard with mainly 
twofold degenerated resonances the advantage of using 



superconducting cavities is obvious. Due to mechanical 
imperfections the degenerated modes show a very weak 
splitting, but nevertheless one is able to resolve all reso- 
nances. As a consequence, all the important characteris- 
tics like eigenfrequencies and widths could be extracted 
with a very high accuracy j|2],fl3| . A detailed analysis of 
the original spectra yields a total number of 1000 - 1200 
resonances for the 2D-billiards (about 660 resonances for 
the circular billiard), and nearly 1900 resonances for the 
3D-billiard. A detailed comparison with numerical data 
confirms that the measured spectra are almost complete 
Q. These eigenvalue sequences {k\, hi, . . . , k n } (with 
k = 2tt/coJ and Co is the speed of light) form the basis 
of the present test of the mode fluctuation distribution. 
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FIG. 1. Investigated billiard systems. In the upper part on 
the left side a Bunimovich stadium billiard (7 = 1.8) and on 
the right a 3D-Sinai billiard are illustrated. In the lower part 
three billiards of the limacon family with different parame- 
ters A are shown changing the chaoticity from nearly regular 
(A = 0.125) to chaotic (A = 0.3). 



IV. APPLICATION TO THE MEASURED 
EIGENFREQUENCY SPECTRA 

In this section we study the mode fluctuation distri- 
bution discussed in Sec. [n] for the experimentally inves- 
tigated systems presented in Sec. III. 



A. 2D-Bunimovich stadium billiard 

First we want to discuss the results for the 7 = 1.8 
stadium billiard. As proven in ]35| ] the stadium billiard 
is strongly chaotic (i.e. ergodic, mixing, K-system). Er- 
godicity, however, does not prevent a system from having 
a family of marginally stable periodic orbits, as long as 
they are of measure zero in phase space. In the stadium 
billiard such a family is given by the bouncing ball or- 
bits (bbo), which have perpendicular reflections at the 
straight line segments. Their contribution dominates the 
fluctuating part of the spectral staircase function |l9| . 
Therefore the above stated conjecture is in its basic form 
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not applicable (see however the refinement given in |l6| ] 
and the discussion below). After unfolding the spectrum 
and calculating the distribution of the normalized mode 
fluctuation W (x) according to Eq. (||) and @), a density 
P(W) results as shown in the left part of Fig. 2. Ob- 
viously it is not a Gaussian, the distribution is shifted 
to positive W(x) due to the existence of the bbo. In 
fact, in |^5| it is shown using the results of J34| that in 
this case W(x) is a bounded function, such that the limit 
distribution P(W) is not Gaussian. 




FIG. 2. Mode fluctuation distribution P(W) for the 
7 = 1.8 Bunimovich stadium billiard. On the left side the 
distribution with bouncing ball orbits (bbo) is shown, which 
is asymmetric with bias to positive W . On the right side 
the distribution without bbo is displayed, which is in very 
good agreement with the Gaussian (dashed curve) expected 
for chaotic systems (Eq. (|lo|)). 

Since P(W) is expected to show Gaussian behavior 
only if the corresponding classical system is strongly 
chaotic excluding stable or neutrally stable motion, the 
contribution of the bbo has to be taken into account in 
the unfolding procedure |l6||, see Sec. ||. Finally, the pre- 
dicted Gaussian distribution of P(W) is obtained, as can 
be seen in the right part of Fig. 2. To have a measure 
how good the distribution P(W) agrees with a Gaussian 
we use the Kolmogorov-Smirnov test, which gives a sig- 
nificance level computed from the maximum value of the 
absolute difference between the two cumulative distribu- 
tions. 

For the case with bbo one obtains a significance level of 
57.2% and for the case with extracted bbo of 75.2%. The 
influence of the bbo which are visible in the distributions 
in Fig. 2 is also reflected in the value of the significance 
level. 

To determine the normalized mode fluctuations, 

Eq. (^J), we have calculated the fluctuating part iVja uc (a;) 

of the unfolded experimental spectra according to 

Eq. (||). For the variance D(x), Eq. (]?]), we have used the 

ansatz D(x) = \J const. ^fx for bbo-contribution included 

after their extraction, 
5 iq6 



1/(2tt 2 ) ln(x) + const. 



and D(x) = 

see Sec. H Then W{x) was calculated for 10° to 
randomly distributed values of x in the interval [0, £ ma x]> 
with a; max being the upper unfolded eigenenergy. The 
constant term in D(x) was chosen according to Eq. (||) to 



give a unit variance for the obtained distribution P(W). 

These results using the experimentally determined set 
of energy levels confirm our previous analysis using the 
numerically computed eigenspectrum u&. 



B. 3D-Sinai billiard 

Another investigated system where the bouncing ball 
orbits play an important role is the 3D-Sinai billiard, 
which is also proven to be strongly chaotic ]39| . In 
contrast to the 2D-Bunimovich stadium billiard the 3D- 
Sinai billiard possesses families of bbo of dimension two 
and three. In Fig. 3 the mode fluctuation distributions 
with and without bbo are displayed. In the distribution 
with bbo (left part of Fig. 3) only slight deviations from 
the Gaussian occur, in particular a significant peak at 
W w — 3 appears. Thus the influence of the bbo is less 
visible in the considered energy range as in the case of 
the distribution for the 7 = 1.8 stadium billiard which 
might be due to the superpositions of different bbo. This 
confirms our previous results of p2^|. Taking into account 
the contribution of the bbo modes one obtains the right 
part of Fig. 3. Calculating the significance level for the 
first case (with bbo) one obtains a value of 78.5% and for 
the second case (with extracted bbo) of 75.2%. 




FIG. 3. Mode fluctuation distribution P(W) for the 
3D-Sinai billiard. On the left side the distribution with bounc- 
ing ball orbits (bbo) is shown, where at W f» —3 a significant 
deviation from the Gaussian occurs. On the right side the dis- 
tribution without bbo is displayed, which shows the predicted 
Gaussian behavior. Again the dashed curve is the expected 
Gaussian for chaotic systems (Eq. (px|)). 



C. Mixed 2D-limacon billiards 

In this section we test the conjecture stated in Sec. || 
for the billiards of the limacon family. We have inves- 
tigated three billiards of different chaoticity with para- 
meters A = 0.125, A = 0.15 and A = 0.3. Investigations 
of the classical Poincare surface of section for these con- 
figurations have shown ]4rj-|4^] that the fraction of the 
chaotic phase space is 55% (A = 0.125), 66% (A = 0.15) 
and nearly 100% (A = 0.3). The quantum mechanical 
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counterpart of these three billiards exhibit the same be- 
havior concerning their chaoticity [|49|J5C| 1 , see also 0] . 

Therefore they are very suitable to study the conjec- 
ture. Note that the classical dynamics of the A = 0.3 
billiard is not completely ergodic, since small stability 
islands still exist plj , see also |52| for analytical results 
around A = 0.25. However the corresponding islands of 
stability are so small that they do not affect the energy 
spectrum in the range under consideration. 




FIG. 4. Mode fluctuations distributions P(W) for the three 
investigated billiards of the limagon family in the upper 
part. All three distributions follow the predicted Gaussian 
for chaotic systems which is represented by the dashed curve 
(Eq. (|lo|)). In the lower part the differences between the his- 
togram and a Gaussian is plotted. 

As can be seen in Fig. 4 the billiard with A = 0.3 shows 
indeed the predicted Gaussian distribution P(W). For 
the calculation of P(W) we followed the same procedure 
described for the stadium billiard with an ansatz for the 
variance D(x) — l/(2ir 2 ) ln(x) + const. . On the other 
hand the two billiards with A = 0.125 and A = 0.15 be- 
long to the class of mixed systems. Their classical coun- 
terparts possess aside from isolated and unstable periodic 
orbits also stable orbits |37|j4^] . From this one would ex- 
pect, that according to this the distribution P{W) should 
show non-Gaussian behavior. However, the histograms 
in Fig. 4, obtained with assuming D(x) cx \fx, which 
gives a good description for D(x) in the considered enery 
interval, allow no significant distinction between P(W) 
and the Gaussian distribution, see also the lower part 
of Fig. 4, where the difference AP between a Gaussian 
and the calculated distribution P(W) is shown. These 
characteristics are also expressed in the significance lev- 
els obtained from the Kolmogorov-Smirnov test which 
lie around 85%. Presumably one needs a large number 
of energy levels to be able to see significant deviations. 



D. 2D-Circular billiard 

Finally we have studied the mode fluctuation distri- 
bution for the circular billiard. Since this system is 
integrable, one would expect a deviation of the mode 
fluctuation distribution from a Gaussian due to the ex- 
istence of neutrally stable periodic orbits. To obtain the 
distribution P(W) we use the ansatz for the variance 
D(x) — const. ■ y/x. In Fig. 5 the distribution P{W) is 
shown, which clearly deviates from the Gaussian. For the 
circular billiard one obtains a significance level of 67.1%. 

This result is in contrast to the result of pi] , where 
no difference between the mode fluctuation distribution 
for the numerical obtained eigcnmodes of a regular sys- 
tem and a Gaussian for chaotic systems could be found, 
although nearly the same number of eigenfrequencies 
(w 660) are used. 
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FIG. 5. Mode fluctuation distribution P(W) for the circu- 
lar billiard on the left side. The histogram shows the expected 
deviation from a Gaussian for chaotic systems which is rep- 
resented by the dashed curve (Eq. (^o|)). On the right side 
one can see the difference between the Gaussian and the his- 
togram. 



V. CONCLUSION 

In this paper we have studied the distribution of 
the normalized mode fluctuations for experimentally ob- 
tained high resolution eigenvalue spectra of several two- 
and three-dimensional superconducting microwave bil- 
liards. The conjecture for this statistical measure states 
that the distribution P(W) of the normalized mode fluc- 
tuations for a given eigenvalue spectrum leads to a Gaus- 
sian when the corresponding classical system is strongly 
chaotic having only unstable and isolated periodic orbits. 

This has been successfully tested using the eigenval- 
ues of a 7 = 1.8 stadium billiard and a 3D-Sinai billiard. 
Both systems are strongly chaotic, but possess families 
of bouncing ball orbits, whose contribution has been sub- 
tracted for the determination of P(W). The same result 
has been obtained with a A = 0.3 billiard of the limagon 
family. Two other investigated billiards of the limagon 
family, the A = 0.125 and A = 0.15 billiard, belong to the 
class of mixed systems. The fact that no visible difference 
from a Gaussian occurs in the distribution P{W) might 
be due to the finite number of eigenvalues in the given 
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energy range we used to calculate the distribution and 
the chaoticities of the systems (55% resp. 66%) which 
are closer to chaotic than to regular dynamics. A varia- 
tion of the included number of modes (500, 600, 1100) 
in the mode fluctuation distribution shows no significant 
change in the results. In the case of the circular billiard 
we get a clear deviation from the Gaussian, as one would 
expect. This is in contrast to the result of 31 1, where no 
such difference for a regular system could be found. 

Therefore characterizing the chaoticity using a small 
number of energy levels with the help of the conjec- 
ture stated in [[L2| is very difficult. For regular resp. 
chaotic systems the reachable number of eigenvalues ob- 
tained by experiments explained in Sec. is sufficient 
to have satisfactory results for the mode fluctuation dis- 
tribution. But for special regular (rectangular billiards 
pl[) or mixed systems the needed number of eigenvalues, 
approximately 10 4 — 10 6 , can experimentally only hardly 
be achieved, so that a practical usage of this conjecture to 
obtain informations about the chaoticity of such systems 
is limited. 
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